Magneto-electronic transport theory in ferromagnets above the Curie temperature 

and in semiconductors 
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Quantitative differences of Lagrange multipliers between standard Fermi-Dirac statistics (FDS) 
and Ionization energy {Ei) based FDS (iFDS) are analyzed in detail to obtain reasonably accu- 
rate interpretations without violating the standard FDS. The resistivity and Hall-resistance mod- 
els in ID, 2D and 3D are also derived to illustrate the transport phenomena in semiconducting 
manganites. It is shown via calculation that the charge carriers in these materials seem to be 
strongly correlated in term of electron-ion attraction or simply, fermions in those materials are 
somewhat gapped due to Coulomb attraction. This Coulomb attraction naturally captures the 
polaronic effect in manganites. Ei is found to be the only essential parameter that predicts 
p{T, doping, pressure, magnetic field) quite accurately. However, this model as will be pointed out, 
is not suitable for metals with free-electrons and strong electron-phonon scattering. It is to be noted 
that iFDS and p(T, doping, pressure, magnetic field) models are only valid in the paramagnetic re- 
gion of ferromagnetic manganites and other doped-semiconductors. Recent X-ray photoemission 
spectroscopy (XPS) studies have indicated that there exists a critical crossover from Mn"^''" to Mn^^ 
depending upon Mn's concentrations in (Gai_a;Mna;)As diluted magnetic semiconductors (DMS). 
Such phenomenon occurring at certain critical concentration directly point towards the applicability 
of Ei based Fermi-liquid theory (iFLT). As such, iFDS is also discussed with respect to DMS above 
the Curie temperature. 

PACS numbers: 71.10.Ay, 75.50.Pp, 72.60.+g, 75.30.Hx 

Keywords: Ferromagnets, Fermi-Dirac statistics, Ionization energy, Resistivity models 



1. INTRODUCTION 



Doped-conipounds including oxides and its elec- 
trical and magnetic measurements have contributed 
enormously to the understanding of electrical prop- 
erties of ferromagnets at T > Tc (paramagnetic 
^ ferromagnetic transition T) and semiconductors. 
The complete mechanism above Tc for ferromag- 
nets is somewhat vague since the variation of 
p{T, doping, pressure, magnetic field (H)) in term of 
hopping activation energy, Ep is still unclear. I.e., the 
variation of Ep with doping is not explicitly predictable. 
Hence, it is essential to study and understand the varia- 
tion of p{T, doping, pressure, H) in order to enhance the 
predictability of electrical properties that may acceler- 
ate the possible applications of these materials. Two- 
dimensional (2D) El based Fermi liquid model was orig- 
inally used to describe c-axis and a6-planes conduc- 
tion of overdoped cuprate superconductors [E S S I3- 
Subsequently, it was further developed to capture both 

^crossover 

(c-axis pseudogap) and T* (spin gap char- 
acteristic T) in a6-planes peculiar conduction involv- 
ing spinous and holons, which is known as the hybrid 
model 0,0. Take note that there are three typographi- 
cal errors in Ref. two of them can be found in Eq. (9) 
of Ref. where {mlml)-^'^ and (27rftVfcB)-3/2 should 
be replaced with (m^m^)^^/"* and (27rfi^/fcB)'^/^ respec- 
tively. Therefore, the fitting parameter, A2 is actually 
equals to {A2/2e'^){mlml)-^/^{2TTh^ /kB)^'"^ ■ One can 



safely claim that manganites are one of the well stud- 
ied compound after Cuprates due to its vast applica- 
tions with hardly any temperature constraints. Literally, 
La-Nd based manganites have been investigated in al- 
most all aspects of experimental techniques, yet there are 
still several theoretical rendezvous on magneto-electronic 
properties both below and above Tc- A wide variety 
of these properties based on doping and Mn's valence 
state in manganites were reported to understand the 
transport mechanism(s) 

Among them, the influence of grain boundary as a bar- 
rier lial. as a region of depleted Tc 0| and polaronic 
effect [ItIIi^ on electrical properties were reported. Di- 
rect proportionality of H with Tc [1^ |2^ |2l| and resis- 
tivity with defects or substrate-film lattice incompatibil- 
ity [Hlilil HI miHlli HI 113 are also regarded as 
equally important to determine the electrical properties 
of manganites. Furthermore, metallic conduction below 
Tc has been studied using double exchange mechanism 
(DEM) between s and d orbitals jsj and the displace- 
ment of hysteresis loop in field-cooled sample with an 
additional scenario of non-linear spin and charge fluc- 
tuations due to magnon |32l |. Explanations in term of 
hopping electrons and DEM [s^, and the influence of 
microstructural transition arises from ionic radius or va- 
lence state of Nd in Nda;Smi_a;Cao.8Mn03 [s^ were also 
reported extensively. 

Presently, the discovery of ferromagnetism (FM) in di- 
luted magnetic-semiconductors (DMS) such as Mn doped 
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GaAs is expected to have an enormous impact on sci- 
entific and technological developments since both elec- 
tronic and magnetic properties are realized simultane- 
ously. Two primary advantages of DMS over mangan- 
ites are the ease of controlling electrons' (n) and holes' 
(p) concentrations and electromigration hardly plays any 
significant role due to zilch oxygen content. Electromi- 
gration in oxides is a process of migration and diffusion 
of mainly oxygen atoms (weakly bound) in the vicinity 
of high temperature, T (random direction), electric, E 
and magnetic fields, H (directional) as well as in high 
vacuum. The attention for the discovery of DMS by 
Munekata et al. [sEl and later by Ohno et al. 
in the early 90s was actually due to the possibility of 
having a brand new type of magneto-electronic materi- 
als other than manganites both in term of experimental 
facts and theories though it does not mimic the excite- 
ment for the discovery of cuprate superconductors in the 
late 80s. Nevertheless, DMS is indeed unique in a sense 
that it exposes the potential applications in a wide va- 
riety of electronic industries namely, high-density non- 
volatile magnetic memory integrated with semiconduc- 
tor integrated-circuits, magnetic sensors, optical isolators 
integrated with semiconductor lasers for optical commu- 
nication systems 38] and nanotechnology [sfll 0| . The 
intensive progress of DMS's sample growth is apparently 
fuelled by the advancement of Molecular-Beam-Epitaxy's 
(MBE) technique in the 90s [illlil lil liil . Basically, 
FM arising from the said doped material was argued 
against the well known Zener DEM 0,0] since DMS re- 
quires long-range spin interaction. This long-range spin 
interaction is necessary to produce FM due to very low 
concentration of magnetic ion (Mn) in DMS and conse- 
quently it is named as such. The highest Curie temper- 
ature {Tc) has been reported to be at 110 K for MBE 
grown, Gao.947Mno.o53As compound [l^l- Recently how- 
ever, Teraguchi et al. '47*1 and Saito et al. [isll have high- 
lighted that Gao.94Gdo.06N and Zno.8Cro.2Te have Tcs 
of about 300 K and 400 K respectively. As a matter of 
fact, it is well known that only oxides namely, Pr,La,Nd- 
Ca,Sr-Mn,Fe-0 compounds provide the established phe- 
nomenon called FM with Tc ranging from above and 
far below room temperature (T,.). In fact, the discovery 
of FM in manganites were first reported by Jonker and 
van Santen |49l Isof in early 50s and antiferromagnetism 
(AFM) in Cd:j;Mn3_^04 by Dey 0. This particular ox- 
ide is also captured in detail by Zener DEM 45. 4 6 j with 
slight improvements in term of polaronic effect [52. Is^ l 
due to Jahn- Teller splitting of Mn'^^ ions. However, FM 
with a Tc up to 110 K that was observed in Ga-Mn- 
As magnetic semiconductors have invoked theoret- 
ical and experimental endeavors due to long-range spin 
interaction. 

Several competing models have attempted to describe 
FM in Mn doped GaAs that includes Ruderman-Kittel- 
Kasuya-Yosida (RKKY) model, DEM and antiferromag- 



netic (AFM) superexchange mechanism as well as the 
Mn-holes (Mn-/i) complex. Actually the RKKY ap- 
proach in the presence of the Friedel oscillations can 
be approximated to DEM for DMS The itinerant 

character of magnetic electrons and quantum oscillations 
(Friedel) of the electron's spin polarization around the 
localized spins were established for the theory of mag- 
netic metals. Hence, the resulting competition between 
FM and AFM interactions give rise to spin-glass freez- 
ing 0, In DMS however, the mean distance be- 
tween the carriers are larger than the distance between 
spins. As a consequence, the exchange interaction me- 
diated by the carriers in DMS favors FM rather than 
spin-glass freezing, which is for mag netic metals 
Originally, Van Esch et al. jH^ IsTl l have given an inter- 
pretation for FM in DMS in term of interaction between 
Mn ions and h surrounding it in which, the Mn ions are 
neutral Mn''+ acceptors associated with 3d^ + h config- 
uration. In addition, this neutral Mn^+ acceptors that 
contribute to magnetic properties could be compensated 
by As, where for a higher concentrations of Mn, instead 
of replacing Ga it will form a six-fold coordinated centers 
with As (Mn*^-^") (5^^J. These centers will eventually 
reduce the magnitude of FM in DMS due to the loss of 
spin-spin interaction between Mn(3c?^) and h. Interest- 
ingly, Dietl et al. have successfully incorporated 
both DEM and AFM superexchange mechanism into the 
Mn-/i complex in order to capture the overall mechanism 
involved in III-V DMS. Such incorporations are impor- 
tant since DEM seems to be the only mechanism that 
predicts Tc accurately after taking band structure effect 
into account [s^ . Issj . DMS was assumed to be charge- 
transfer insulators by Dietl et al. because h from Mn^+ 
does not reside on d orbitals but occupies an effective 
mass Bohr orbit with large intrasite correlation energy. 
Moreover, Mn ions also function as a source of localized 
spins as well as acceptors [H^ . Dietl et al. '5?, 'SS^ have 
further suggested that FM between localized spins {S — 
5/2) of Mn is mediated by delocalized or weakly localized 
holes surrounding it and the resulting spin-spin interac- 
tion are indeed long range as it should be. On the other 
hand, the short-range AFM superexchange (if any) arises 
from the spin polarization of occupied electron bands (va- 
lence band). As such, AFM superexchange is mediated 
by short-ranged localized spins unlike FM where it re- 
quires at least weakly localized spins for long range in- 
teraction. The readers are referred to Dietl et al. |54| for 
a thorough review on DMS. 

Apart from that, Okabayashi et al. |H3 have pointed 
out convincingly that there exists an indeterminacy on 
pin pointing whether Mn^+ (negative) or Mn^^^ (neu- 
tral) ion that dominates Gao.93Mno.07As DMS. Both ions 
(Mn^+ and Mn'^+) were found to be favorable towards 
FM. In order to explain this scenario, Ando et al. jH^ 
have proposed that Mn'^+ ions populate at very low Mn 
concentrations due to low screening effect or large long- 
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range Coulomb interaction. I.e. long-range Coulomb in- 
teraction will be screened-out for higher Mn concentra- 
tions. This interesting effect of transition between Mn^+ 
(negative) ^ Mn^+ (neutral) with Mn doping can also 
be analyzed in accordance with iFLT. The mechanism 
that governs DMS above Tq and the predictability of its 
/9(r, doping) and i?/f(r, doping) as well as the changes in 
Tc with H are also evaluated. Actually, the transition of 
Mn^+ (negative) ^ Mn'^+ (neutral) can be described by 
simply calculating Ej for each doping from its respective 
p{T, doping). 

In this paper, iFDS is re-derived to extract the La- 
grange multipliers so as to trigger sufficient interest for 
applications in other compounds such as ferromagnetic 
manganites and DMS above Tc- p{T) curves are sim- 
ulated at various doping or gap {Ei), pressure and T 
to further enhance its applicability. Basically, we will 
stress the polaronic effect in semiconducting mangan- 
ites and DMS in order to obtain accurate predictions 
of the charge carriers transport mechanisms above Tc 
where the changes in Ei is well accounted for with dop- 
ing. In addition, it is also shown with detailed derivation 
of somewhat different Lagrange multipliers that sepa- 
rately influence the standard FDS and iFDS. In addition, 
the probability functions of electrons and holes, charge 
carriers' concentrations, p(r, doping) and _R^f(r, doping) 
models in ID, 2D and 3D are derived as well via the 
standard quantum statistical method. Interpretations of 
electrical properties based on these models for mangan- 
ites and DMS are highlighted in detail. Quantum statis- 
tical method is used instead of other approaches because 
the former method will lead to the final results to un- 
derstand p{T, doping, pressure, H) without imnecessary 
approximations. 



2. THEORETICAL DETAILS 

2.1. Fermi-Dirac distribution function based on 
ionization energy 

Both FDS and iFDS are for the half-integral spin par- 
ticles such as electrons and holes. Its total wave function, 
^ has to be antisymmetric in order to satisfy quantum- 
mechanical symmetry requirement. Under such condi- 
tion, interchange of any 2 particles {A and B) of differ- 
ent states, ipi and ipj {j ^ i) will result in change of sign, 
hence the wave function for Fermions is in the form of 



respectively. Equation |^ is nothing but Pauli's exclu- 
sion principle. The one-particle energies Ei, E2, E3, 
Em for the corresponding one-particle quantum states 
qi, 92, 93, qrn Can be rewritten as {Einuiai state ± 

{Einitial state ± ^1)2, {Einitial state 
{Einitial state ± Ei)m- NotC here that Einitial state + 

El — Eeiectrons and Einitial state ^ Ej = Ehoies- Sub- 
sequently, the latter {Einitial state ± Ei)i version where 
i — I, 2, 3, m with Ej as an additional inclusion 
will be used to derive iFDS and its Lagrange multipliers. 
This ±Ei is inserted carefully to justify that an electron 
to occupy a higher state N from initial state M is more 
probable than from initial state L if condition Ei{M) < 
Ej{L) at certain T is satisfied. As for a hole to occupy a 
lower state M from initial state N is more probable than 
to occupy state L if the same condition above is satisfied. 
Einitial state IS the energy of a particle in a given system 
at a certain initial state and ranges from +00 to for elec- 
trons and to —00 for holes. In contrast, standard FDS 
only requires Ei [i — 1, 2, 3, m) as the energy of a par- 
ticle at a certain state. Denoting n as the total number 
of particles with ni particles with energy {Einitial state ± 
Ei)i, 712 particles with energy {Emtttai state ± Ei)2 and 
so on implies that n = ni -I- 722 + f^3 + ••■ + nm- As a 
consequence, the number of ways for qi quantum states 
to be arranged among ni particles is given as 



P{ni,qi) 



91 ! 



ni!(gi - ni)! 



(2) 



Now it is easy to enumerate the total number of ways 
for q quantum states {q = qi + q2 + q^ + ■■■ + 9™) to 
be arranged among n particles, which is 



P{n,q) = \{ 



^Jl n.i\{qi - ni)\ 



(3) 



The most probable configuration at certain T can be 
obtained by maximising P{n, q) subject to the restrictive 
conditions 



i 

00 

drii = 0. 



(4) 



*.,j(Ca,Cb) = i>^{CA)^,{CB) ~ ^.(Cs)^,(C^), (1) 

The negative sign in Eq. 0J that fulfils antisymmet- 
ric requirement is actually due to one of the eigenvalue 
of exchange operator |20|, P = — 1. The other eigen- 
value, P = -1-1 is for Bosons. Ca and Cb denote all the 
necessary cartesian coordinates of the particles A and B 



i 

00 



initial state 



initial state 



± Ei\n, = E, 



± Ei\dn^ = 0. 



(5) 



The method of Lagrange multipliers |60j can be em- 
ployed to maximise Eq. Hence, a new function. 
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F{xi,X2,X3, ■■■fi, A, ...) = / + Ai/i + A/2 +... is introduced 
and all its derivatives are set to zero 



^=0; ^^0: ^=0. 



(6) 



As such, one can let the new function in the form of 



ial state 



After applying Stirling's approximation, dF/drii can 
be written as 



V 



V 





■-A?l2k2" 


j exp 


2m 



dk, (12) 





■-a;i2j^2- 


/k^exp 


2m 



dk. 



(13) 



n is the concentration of electrons whereas p represents 
holes' concentration. V denotes volume in k-space. The 
respective solutions of Eqs. l(T^ and ((T^ are given below 



dF 

drii 



\n{qi - Tii) - Inrii + /i + X{Ei„iUai state ± Ei)i 



= 0. 



(8) 



Thus, the Fermi-Dirac statistics based on ionization 
energy is simply given by 



(9) 



qi exp[^ + \{E^nitial state ± Ei)i] + 1 ' 



The importance of -E/'s inclusion is that it can be in- 
terpreted as a charge gap that will be described later and 
also, particularly the Ej can be used to estimate the re- 
sistivity transition upon substitution of different valence 
state ions. 



2.2. Lagrange multipliers 

By utilizing Eq. I^l and taking exp[^ + \{E ± Ej)] > 
1, one can arrive at theprobability function for electrons 
in an explicit form as |3] 



/e(k) = exp 



-,-X[^+E, 



(10) 



Similarly, the probability function for h is given by 



exp 



2m 



(11) 



The parameters n and A are the Lagrange multipliers. 
H = h/2TT, h = Planck constant and m is the charge 
carriers' mass. Note that E has been substituted with 
S2k2/2m. In the standard FDS, Eqs. ^ and ^ are 
simply given by, /e(k) = exp[-^-A(?i2k2/2m)] and /^(k) 
= exp[/i -|- A(fi,2k2/2m)]. Equation can be rewritten 
by employing the 3D density of states' (DOS) derivative, 
dn = V'k'^d'k/2Tr'^, that eventually gives 0| 



Me 



fJ-h 



In 



In 



n /27rAfi2 
V V m 



V V m 



3/2' 



3/2' 



(14) 



(15) 



The subscripts e and h represent electrons and holes 
respectively. Separately, Eq. (O after neglecting ±Ei 
(for FDS) can be written as 



E 



— e ^ 
2A 



k'' exp 



3/2 



2m 



27rAft2 



(16) 



Quantitative comparison between Eq. H16(l and 
with the energy of a 3D ideal gas, E = 3nkBT/2, af- 
ter substituting Eq. H14|) into Eq. (|16|) will enable one to 
conclude Xfds — ^/ksT. fcs is the Boltzmann constant. 
Applying the identical procedure to iFDS, i.e. by em- 
ploying Eqs. Unj) and (jlTJ, then Eqs. and (O are 
respectively rewritten as 







V 



P = 2^' 



2^-^^^ J k2exp(^A^^dk. (18) 

— oo 

The respective solutions of Eqs. H17(l and 1)18(1 are 



fi + XEj = -In 



V \ m 



(19) 
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fi- XEi 



In 



3/2' 



(20) 



Note that Eqs. (f11?|l and simply imply that 

lie{iFDS) ^ lie + XEi and iih{iFDS) ^ Hh - XEj. 
Furthermore, using Eq. one can rewrite Eq. H16(l as 



P = 



f,{E)N,(E)dE, 



fh{E)Nh{E)dE. 



(24) 



(25) 



E 



3V 



k exp —A 



2m 



2A V27rA;i2y 



3/2 



(21) 



Again, quantitative comparison between Eq. (|21|l and 
with the energy of a 3D ideal gas, E = 3nkBT/2, after 
substituting Eq. 119|l into Eq. H21|l will enable one to 
determine A. It is found that A remains the same as 
l/ZcsT. I.e., Xfds = Xipos as required by the standard 
FDS. Hence, the relationship between FDS and iFDS in 
term of Lagrange multipliers has been derived and shown 
clearly. 



Existence of Eg , which is the energy gap due to energy 
band splitting or lattice based gap is not inserted explic- 
itly thus it is (if any) can be coupled with Ej, which is 
tied to ions via Coulomb attraction. Having said that, 
now it is possible to obtain the geometric-mean concen- 
trations of electrons and holes for ID, 2D and 3D respec- 
tively in the forms of (assuming n « p) 



np(lD) = 



2-K 



X exp 



-Ei 



1/2 



(26) 



2.3. Resistivity models 

Denoting /i = —Ep (Fermi level), A — X/ksT, 
h'^k'^/2'm = E and substituting these into Eqs. ((TUIl 
and pi|l will lead one to write 



feiE) ^ exp 



Ep — Ej 



E 



(22) 



1^{2D) 



h) exp 



7^^(31?) = 2 



In which 



knT 



2Trff 
X exp 



3/2 



-Er 



knT 



(m>;;)3/4 



'Ei 



knT 



(27) 



(28) 



fhiE) = exp 



'E-Ei-Ep' 
k^ 



(23) 



At this point, one might again wonder the reason for 
-B/'s inclusion. The unique reason is that it directly de- 
termines the kinetic energies of electrons which carry the 
identity of its origin atom. Detailed experimental impli- 
cations are given in the discussion. These iFDS probabil- 
ity functions for electrons and holes are unique in a sense 
that it allow the prediction of charge carriers' concentra- 
tions at various T and doping. It is worth noting that, 
—Ei in Eq. H23(l for h follows naturall y fr om the Dirac's 
theory of antiparticle interpretations [61|. Besides, the 
charge carriers are not entirely free since there exist a 
gap-like parameter that can be related to electrons-ion 
or Coulomb attraction. In fact, application of Eqs. H22|l 
and H23|l in semiconducting manganites will be explained 
in section IIIB. The general equations to compute charge 
carriers' concentrations are stated below. 



71(31?) 



'keT' 

2TTh^ 



3/2 



3/2 



X exp 



Ep — Ej 



p{3D) 



kuT 



knT 



-1 3/2 



(29) 



2T:h? 



X exp 



—Ep — El 



(30) 



It is apparent that n(lD,2D,3D) cx exp[(i?i? — 
Ei)/kBT] and p(lD,2D,3D) cx exp[(-SF - Ei)/kBT] 
while the rest are just typical physical constants ex- 
cept T. As such, the constant, Ep in the exponential 
term will appear where appropriate if one assumes n 3> 
p or p > n instead of ^np. The DOS, N{E, ID) = 
[E-^l'^yf^^^)l-Kh, N{E,2D) = m*/nh'^ and N{E,3D) 
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= (£'^/^/27r^)(2m*//i^)^/^ were employed in which, m* is 
the effective mass. Consequently, the resistivity models 
for ID, 2D and 3D can be derived from p = m/ne^T by 
taking l/r = ^£1=1,2,3^^ as a consequence of electron- 
electron scattering and m = ml « mjjj ~ (m*mj^)^/^. 
The subscript D = 1,2,3 represents dimensionality, i.e. 
j4i, A2 and A3 are T- independent scattering rate con- 
stants in ID, 2D and 3D respectively, r denotes scatter- 
ing rate due to electron-electron scattering in the absence 
of H. Therefore, the respective p{T,Ei) are given by 



p{lD) = 



1/2 



exp 



Ei 



(31) 



The subscripts x, y and z represent the axes in x, y and 
z directions while the scattering rate, 1/th = ^c22 3-^^ 
in which A]jJ^2 3 niay not be necessarily equals to AD^2.■i^ 
though both Ad^2,3 and J2 3 independent of T. 
A'--^^ and th denote the T-independent scattering rate 
constant and scattering rate respectively with applied H. 
In a steady state of a static H and E, dvz/dt ~ dvy/dt 
= and Vz — hence, E, can be obtained from 



E, 



eHx^yTH 



(38) 



In addition, it is further assumed that Px{T) — py{T) — 
p,[T) = p{T). is defined as EjjyH,, jy = Ej,/p(r) 
and tan^^^■' = E^/E^. Parallel to this. 



p{2D) 



A^nh^ 



Texp 



Er 



(32) 



(,) tan4"V(r) 

- 7? • 



(39) 



3/2 



xT^/^exp 



Ei 



knT 



In which 



ks 

xT^/^exp 



Phi'SD) 



A3 \2T:h^ 
2? 

xT^/^cxp 



3/2 

Ej — Ep 
keT 

3/2 

Ei + Ef' 



knT 



2.4. Hall resistance 



(33) 



(34) 



(35) 



The equations of motion for charge carriers under the 
influence of static H and electric field (E) can be written 
in an identical fashion as given in Ref. ~62\, which are 
given by 



m 



d 
dt 

d 
di 



I 

TH 

1 

TH 



Vz = -eE^ + eHxVy. 



(36) 



(37) 



jy is the current due to charge carriers' motion along 
y-axis and fl^^* is the Hall angle. Furthermore, tanfl^-* 
can be rewritten as —eilxTH/m. Therefore, it is easy 
to surmise that cot^^' cx T^. After employing Eq. H32() 
with Ep {n ^ p) and Eq. (|34() . then one can respectively 
arrive at 



R' 



{2D) 
■H 



AoTrh"^ 



A^^ emlkB 



exp 



Er — Ep 



knT 



, (40) 



R 



{3D) 
H 



A. 



2A^ 

,y-3/2 



i(m*)3/2 



2tiW- 



21-3/2 



exp 



E^ 



ks 



knT 



(41) 



The negative charges in Eqs. (|38|l . 140() and H41|l are due 
to the assumption that the charge carriers are electrons. 
Note that Rh{T,1D) for any given samples that exhibit 
purely ID conduction is obviously irrelevant or simply, 
could not be derived with above procedures, since Hall 
effect requires at least 2D conduction. 

2.5. Special cases 

There are nowhere in the above derivations that take 
into account any free electrons and T-dependence of 
electron-phonon scattering. Hence, the models derived 
thus far are obviously not suitable for such applications 
except for semiconducting free electrons above conduc- 
tion band. In this case, Eqs. H24|l and (|25|l should be 
integrated from iJg — s- 00 and — > — 00 respectively after 
replacing i?/ = in Eqs. and As for metals 
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with free electrons and strong phonon contributions, it is 
advisable to switch to the well known Bloch-Griineisen 
formula given by |6^ , 



p{T,ZD) = po + At 



1287rm*(fcsT)5 



sinh^ 



-dx. 



(42) 



\tT — electron-phonon coupling constant, po — p{T = 
0), TO* — average effective mass of the occupied carrier 
states, Qd — Debye temperature, n = free electrons con- 
centrations. As a matter of fact, one should not be en- 
couraged to substitute any of the Eqs. if^ — for n 
into Eq. just to capture the electron-phonon scat- 
tering because the scattering of free electrons considered 
in Eq. H42|l may not be compatible with (7appe(i-electrons' 
scattering of iFDS, unless proven otherwise. 



2.6. The rationale behind iFDS 

Recall that Ei is the gap from electron-ion attraction 
along the easiest path while Eg is the averaged gap from 
the energy band splitting. Eg from the theory of elec- 
tronic band structure (TEBS) does not give the correct 
prediction on how certain ionic substitutions affect the 
transport properties of strongly correlated systems. As 
I have pointed out in the last paragraph of the intro- 
duction, iFDS follows a brand new approach instead of 
the common TEBS, which is not suitable and/or reliable 
for the strongly correlated systems namely, superconduc- 
tors, ferromagnets and ferroelectrics. In fact, application 
of iFDS in the normal state of high- Tc superconductors 
have been carried out while its consequences in ferro- 
electrics has been reported as well Usually, one 
needs to incorporate TEBS into FDS to take the lattice 
properties into account to expose the transport mecha- 
nism and one can also include the defects so as to arrive 
at a comprehensive theory. However, the magnitude of 
lattice effects and anisotropies calculated from TEBS is 
obviously not suitable for the transport theory because 
the charge carriers move along the easiest path comply- 
ing the principle of least action. This is why TEBS-FDS 
is not reliable when it comes to the transport properties 
of strongly correlated systems. It is also important to 
realize that for simple systems, the lattice properties are 
isotropic therefore the magnitude of say. Eg calculated 
from TEBS is comparable with Eg from the easiest path, 
thus simple systems are very well described by TEBS- 
FDS. Unlike TEBS-FDS, iFDS takes care of the crys- 
talline lattice via Ej in accordance with the easiest path 
(transport measurements) since any changes to lattice's 
crystallinity due to substitution will be mirrored through 



the easiest path's Ej and subsequently shall be picked 
up by the electrons, influencing its transport mechanism. 
Whereas, the TEBS-FDS theory assumes that all the lat- 
tice effects including Eg and its magnitudes are isotropic, 
which will be imposed on the electrons unnecessarily. In 
other words, it is logically desirable to let the electrons 
reveal the obstacles via iFDS rather than forcing Eg onto 
the electrons in order to calculate the electronic transport 
properties. 

I would like to stress here that if one substitutes Ei 
with Eg in Eq. H33|l to comply with the standard FDS, 
then one has to feed the value of Eg from elsewhere 
(such as TEBS) in order to predict the transport prop- 
erties and apparently, TEBS is not suitable for the rea- 
sons described above. Hence, the Eqs. H26(l — H35|l and 
Eqs. (|40|) ~ H41|) can be directly employed to analyze the 
doping effect on the transport properties of ferromagnets 
above Tc, be it diluted (DMS) or concentrated (mangan- 
ites). Take note that the proposed theory based on iFDS 
is in its infancy, thus the correct method to insert defects 
via El will be developed later on. Consequently, the ap- 
plication of iFDS is to be confined to pure crystalline and 
polycrystalline samples. The above-stated arguments do 
not invalidate the importance of TEBS since the lattice 
properties and its consequences on electronic structure 
are indeed invaluable to derive lattice based Hamiltonian. 



3. DISCUSSION 
3.1. Simulated curves 

Figures^, lb and Ic illustrate the variation of p{T) 
from Eqs. (|31|) — H33|) respectively with conduction dimen- 
sionality and doping parameter {Ei). The proportion- 
ality of Eqs. (|3U-(|33J with Ei is similar to Eqs. 
and 1)35(1 . thus it does not matter whether one chooses 
Y^np, n or p in the subsequent analysis as well as for man- 
ganites. One can also identify the p{T) transition from 
metallic semiconducting conduction with increasing 
Ej. It is also worth noting that p(r,lD), p(r,2D) and 
/9(T,3D) are cx T^/^, T and y/T respectively if and only 
if El <^ T. Another point worth to extract from these 
curves are the relationship between Tcrossover and Ej, 

where Tcrossover < Ej for ID, Tcrossover = Ej for 2D 

and Tcrossover > Ej for 3D. Apparently, these relations 
are again due to the proportionalities of p(T,lD) cx T'^/^, 
p(r,2D) cx T and p(T,3D) cx %/T. Figures Hi and 2b 
plot the simulated Rh{T) curves in 2D and 3D as well 
as at different Ei{0 K, 150 K, 310 K) that follow from 
Eqs. (|40|l and (|41|l . There are no significant differences of 
Rh between 2D and 3D since i?//(2D) and i?H(3D) are 
cx to T~^ and r~3/2 respectively. Besides, the T from 
exp[i?//r] also inversely proportional to both i?/f(2D) 
and (3D). These scenarios will always lead Rh to in- 
crease with lowering T without any observable Tcrossover 
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regardless of Ej's magnitude, unlike p{T). It is conve- 
nient to directly quantify p{T) variation with doping by 
relating Ej as a doping parameter, as will be discussed 
in the following paragraph for manganites at T > Tq. 



3.2. Manganites 

As stated previously, Jonker and van Santen [491 ISC 
suggested that FM is due to indirect coupling of d-shells 
via conducting electrons. Subsequently, Zener 0, 0| , 
Anderson and Hasegawa |^ have provided sufficient 
theoretical backgrounds on DEM. However, this paper 
will not discuss the property of DEM below Tc, instead 
the electrical properties above Tc (paramagnetic phase) 
will be addressed in detail in which, DEM is employed 
at T < Tc (ferromagnetic phase). It is interesting to 
observe reduced p{T) and increased Curie temperature 
(Tn) at higher H for Lai_2;Ca2;-Sra;Mn03 compounds 
[Ha. IstI IbsI . The results that larger H giving rise to con- 
ductivity at r > Tc point towards the enhancement of 
conductivity from DEM where the exponential increase 
of p{T) is suppressed with larger H. Additionally, small 
polarons have also been attributed [g^ to play a signifi- 
cant role on p{T) at T > Tc- As a matter of fact, this po- 
laronic effect is naturally captured by Eqs. 1)33(1 — 1)35(1 in 
which the gap-parameter, Ej that represents electron- ion 
attraction is also a parameter that measures the combi- 
nation of electrons and its strain field due to neighboring 
ions, which is nothing but polarons. The absolute value 
of Ej can be obtained from 



Er 



(43) 



e and eo are the dielectric constant and permittivity 
of free space respectively, rs is the Bohr radius. Fur- 
thermore, the decrement of Ej with H indicates that rs 
increases with H. Identical relationship was also given be- 
tween polaronic radius, and Ep by Banerjee et al. [70j . 
Actually, Millis et al. [53j have somewhat proved the in- 
adequacy of DEM alone to describe p(T, doping,!!) and 
reinforced the need to include small polarons as a conse- 
quence of Jahn- Teller (JT) splitting of Mn^^^ ions. This 
statement was furtherjustified by experimental work of 
Banerjee et al. [t^, EI Iz3 in which, they have estab- 
lished the existence of small polarons at T > Tc in 
Lao.5Pbo.5Mni_a;Cr2,03 for a; = ^ 0.45 using ther- 
moelectric power and positron annihilation lifetime mea- 
surements. Banerjee et al. suggested that the substi- 
tution of Cr'^^ into Mn sites localizes electrons that 
gives rise to p(T) III [73. However, the increment 
of activation energy, Ep with x is still unclear in term 
of Cg electrons' localization. On the other hand, p{T) 
(Eqs. (|2Sll-(IS3) based on iFDS could explain the in- 
crement of Ej with X as well as the structural changes 



accompanied by this doping, which is due to the fact that 
valence state of Cr and Mn may change with doping that 
can be calculated with Eg. I44II as will be shown later. 

In addition, Moskvin |7J,l75j reinforces the importance 
of considering different charge distribution in MnOe, Mn 
and Oxygen instead of just considering the DEM. Par- 
allel to this, doping- friendly resistivity models derived 
from iFDS will indeed lead to identical conclusions with 
Billinge et al. and Moskvin since iFDS naturally consid- 
ers polarons and the substitution of Ca^+ or Cr^+'^+'^+ 
into Lai_a;Caa;Mni_yCrj,03 system does indeed modify 
the overall charge distribution in accordance with the va- 
lence state of Ca, La, Mn and Cr. Add to that, Louca 
and Egami |^ invoked the JT distortion to describe the 
effect of lattice on Tc . They have utilized the results of 
pulsed neutron-diffraction experiments to conclude that 
the variation in Mn-0 bond length with Sr substitution in 
Lai_a;Sr2;Mn03 compound can be related to JT. Again, 
note here that the change of Mn-0 length with Sr substi- 
tution implies the valence state of Mn varies with doping. 
In an identical compound of Lai_a;Caa;Mn03, substitu- 
tion of Ca into La, will have to satisfy the inequality of 
average Ej between Ca2+ (867 kJmol'^) and LaP+ (1152 
kJmol^^) i.e., i?/(La'^+) > i?/(Ca^+). Here, one can eas- 
ily fix the valence state of Ca^"'' and La^^^ as noted. As 
a consequence, this will ease the prediction of p(T) with 
doping. If one of the ions is multivalence, then the linear 
algebraic equation as given below must be used to pre- 
dict the valence state of the multivalence ion from p(T) 
curves 



S 1 ^ 1 

- V En + -Y^En = -Y^En 

■' i=z+l i-1 ^ "-I 



(44) 



Eji above has i — z + I, z 



z+l 

2,..., z + j and j — 1, 2, 3,.... It is solely due to 



The first term, i 

id j = 

multivalence ion for example, assume Mn^+'^+ is substi- 
tuted with Nd'^"'" (Lao.7Cao.3Mni_j;Nda;03) hence from 
Eq. ((44(1 . the first term is due to Mn'*+ ion's contribu- 
tion or caused by reaction of the form Mn'^^ — electron 

Mn4+ (4940 kJmol'^), hence 2 is equals to 1 in this 
case and 8 represents the additional contribution from 
Mn^'^. The second (i = 1, 2, 3, z) and last (i = 1, 
2, 3, q) terms respectively are due to reaction of the 
form Mn — 3(electrons) — > Mn'^+ and Nd — 3(electrons) 

Nd3+. Recall that g = z = 3-h and i = 1, 2, 3,... 
represent the first, second, third, ... ionization energies 
while j — 1, 2, 3, ... represent the fourth, fifth, sixth, 
... ionization energies. Therefore, z -\- 5 gives the min- 
imum valence number for Mn which can be calculated 
from Eq. (|^ . Now, it is possible to explain the doping 
effect in Lao.5Pbo.5Mni_a;Crj;03 system 

IzSEIl- The 

inequalities of Ejs are given as Mn'^^ (1825 kJmol"^) 
> Cr3+ (1743 kJmol-i) and Mii^+ (2604 kJmol"!) > 
Cr''+ (2493 kJmol^^). These relations strongly indicate 
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that p{T) should decrease with Cr^+ content contradict- 
ing with experimental data from Refs. [toI FH . The only 
way to handle this situation is to use Eq. H44|l so as to 
calculate the minimum valence state of Cr'^+'^, which is 
3.033+. Actually, the valence state of Cr that substitutes 
Mn'^+ is Cr>'^ °'^^ and of course, the valence state of Mn 
is fixed to be 3+. There is no need to var y it because 
p{T) was found to increase with Cr content jz^^Sl- Ac- 
tually, the existence of Mn^+ is unlikely and negligible, 
while Mn^"*" (if any) will further increase the value of 5 
in Cr'^+*, in which both scenarios comply with iFDS via 
Eq. and the resistivity measurements [t^, 0| . 

The effect of hydrostatic (external) pressure (P) and 
chemical doping (internal P) on metal-insulator tran- 
sition of Pro.7Cao.2Sro.iMn03, Pro.7Cao.2iSro.o9Mn03, 
Pro.58Lao.i2Cao.3Mn03 and Pro ^iLa,n ifiCan 3Mn03 sys- 
tems were reported by Medvedeva et al. [22|. It is found 
that p(T) and Tc are observed to be decreased and in- 
creased respectively with increasing P ranging from — > 
15 kbar. As anticipated, p{T) above and below Tc have 
been decreased significantly with P i.e. P affects both 
Mn-O-Mn bond angle and length. Hence, it is apparent 
that doping and P give rise to the variation in the valence 
state of Pr and Mn in order to achieve a certain crys- 
tal structure and simultaneously increase the number of 
charge carriers. As such, changes in p(T) above Tq can be 
very well accounted for with Eqs. (I2Hll-(|Sni), 
where P reduces Ej (increases tb) of certain ions in a 
similar fashion to doping (internal P). However, this pa- 
per as stated earlier does not attempt to describe the cor- 
rect mechanism(s) involved below Tc and its variations 
with P and doping. In short, at T < Tc, the interac- 
tions among DEM with JT or polarons or all rnay play 
a significant role as suggested by Billinge et al. and 
Medvedeva et al. [t^ . 

3.3. Diluted magnetic semiconductors 

It is quite straight forward to extend iFLT described 
above to DMS. Firstly, Oiwa et al. Matsukura et 
al. jz^i Dietl et al. [S0| and lye et al. have done 
detailed investigations on Ga-As-Mn based DMS. The 
variation of electrical resistivities with temperature and 
doping for Gai-^jMn^jAs in Refs [ll [ll H^l (sample 
Al) differ slightly from Ref [HI (sample B2). It has 
been reported that the maximum Tc is ~ 70 K for 
Gao.957Mno.o43As |i3 whereas Gao.947Mno.o53As had a 
Tc of ^ 100 K. Apart from Tcs, there is a switch-over 
in the variation of p{T, x) where p(T, x) decreases with 
X initially before switching to increasing p(T, x) with 
X. p(r, x)'s variation for sample Al is p(r, 0.015) > 
p(r, 0.022) > p(T, 0.071) > p(T, 0.053) > p(r, 0.043) > 
p{T, 0.035) while for sample B2, p{T, 0.015) > p(T, 0.022) 
> p(r, 0.071) > p{T, 0.035) > p(r, 0.043) > p(r, 0.053). 
The initial reduction of p(T) with x is due to Ej (Mn'^+ 



= 1825 kJmol-^) < Ei {Ga.^+ = 1840 kJmol"!), i.e. it 
supports that Mn^+ substitutes at Ga'^+ sites. How- 
ever, p(T)'s switch-over (from decreasing to increasing 
with x) at critical concentrations, Xci = 0.035 and Xc2 
— 0.053 for samples Al and B2 respectively can only 
be explained if one employs the mechanism proposed by 
Van Esch et al. |5^ [53 and Ando et al. [53 in which, 
y[j^6As formations is substantial above Xc in such a way 
that Mn'^+ ions do not substitute Ga'^+ ions. Therefore, 
p{T) will be influenced by Mn^^^ clusters, defects and 
Ga-Mn-As phase simultaneously significantly above Xc. 
This is indeed in fact in accordance with iFDS based re- 
sistivity models since if one assumes Mn^"'' substitutes 
Ga^+, then p(T) should further decrease with x with- 
out any switch-overs. Hence, the formation of Mn^'^* 
with X is inevitable in III-V Gai_a;Mnj;As DMS. Again, 
Tc of DMS is also directly proportional to H, identical 
to manganites. All the values of Ej discussed above (in- 
cluding in section IIIB) were averaged in accordance with 
i?/[Ar^+] = X]i=i "T^ ^^"^ should not be taken literally 
since those Ejs are not absolute values. The absolute 
values need to be obtained from the rs dependent Ej 
equation stated earlier [sj [s^ . Prior to averaging, the 
V*, 2"'', 3'"'* and A*^ ionization energies for all the el- 
ements mentioned above were taken from Ref. 84] and 
the predictions stated above are only valid for reasonably 
pure materials without any significant impurity phases as 
well as with minimal grain boundary effects. 

4. CONCLUSIONS 

In conclusion, the ionization energy based Fermi-Dirac 
statistics is useful to estimate the transitional progress of 
p{T, doping, pressure,!!) from metallic — > semiconduct- 
ing or vice versa in both manganites and diluted magnetic 
semiconductors. This is made possible by an additional 
unique constraint, which is nothing but the ionization en- 
ergy that captures the electron's kinetic energy and maps 
it to its origin atom. The relation of Lagrange multipli- 
ers (A and /i) between FDS and iFDS have been derived 
explicitly solely to fiush-out any misinterpretations that 
will lead to further complications in describing experi- 
mental data on semiconducting ferromagnets and DMS 
at T > Tc. The presented iFDS model however, does 
not admit completely free-electrons and strong electron- 
phonon scattering. Importantly, Ej captures the pola- 
ronic effect quite naturally to explain the paramagnetic 
phase electrical properties of the said materials accu- 
rately. 
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FIG. 1: Simulated p{T,Ei) curves from Eqs. (glj, JHU 
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respectively. All curves for both dimensionalities (2D and 
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K. Note that Rh{T, Ej) curves are inversely proportional to 
T regardless of the magnitude of Tcrossover that can be easily 
verified from Eqs. 14011 and 14111 . Therefore, one could not 
excerpt sufficient quantitative differences between 2D and 3D 
Rh{T,Ei) as opposed to p{T,Ei). 
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